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摘 要 ASAA E (LBE) 是 一 种 性 能 优异 的 材料 ， 在 加 速 驱动 系统 (ADS) 中 作为 散 裂 攻 和 冷却 剂 材 料 . 然而 , 在 
ADS 中 ,贮存 LBE 的 容器 和 结构 钢材 料 在 高 温 下 与 LBE 直接 接触 时 会 被 严重 腐蚀 . 本文 研究 该 腐蚀 和 金属 Cr, Ni 减 
轻 该 腐蚀 的 机 制 . 研究 发 现 ， 铁 在 LBE 中 被 腐蚀 的 微观 机 制 是 由 于 铁 原 子 与 LBE 原子 相互 渗透 ,破坏 了 铁 材料 结构 ; 
铁 中 掺 入 Cr, Ni 可 以 有 效 减轻 铁 的 腐蚀 , 其 微观 机 理 为 金属 Cr、 Ni 加 入 到 铁 后 阻止 了 LBE 原子 渗入 铁 材料 中 破坏 其 结构 . 
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Abstract The liquid-bismuth eutectic (LBE) is an excellent material as coolant and spallation target 


material in the accelerator driven system (ADS). However, the molten metal has a severely corrosion 


for structural steel material and the container storing LBE, especially at high temperatures. This 


paper studies the mechanism of the corrosion and mitigating the corrosion by using Cr and Ni. It is 


concluded that the micro-mechanism of the corrosion of the iron in the LBE is due to the mutual 


penetration between the iron and the LBE atoms destroys the structure of the iron. Meanwhile, after 


adding the Cr and Ni atoms in the iron, the corrosion of the iron is effectively reduced, which owes to 


the mixing of the Cr and Ni atoms in the iron prevents the LBE atoms from penetrating into the iron 


material and destroying its structure. 


Key words 
05| F 


LBE 是 ADS 中 一 种 性 能 优异 的 冷却 剂 和 散 裂 
靶 的 材料 由 3， 这 是 因为 其 低 熔 点 (398.25 K), mb 
点 (1943.15 K) 回 ， 低 蒸汽 压力 和 良好 的 传 热 性 能 A. 
这 些 优 异 的 性 能 使 得 ADS 运行 过 程 更 加 可 靠 和 安 
全 D, 此 外 , LBE 不 仅 与 水 和 氧气 成 化 学 惰性 ， 而 
且 每 1 GeV 的 能 量 能 使 其 产生 30 个 中 子 回 ， 与 其 
他 材料 相 比 ， 能 够 产生 更 多 的 中 子 。 基 于 以 上 优异 
性 能 使 得 LBE 成 为 了 ADS 的 最 佳 材 料 之 一 。 

众所周知 , 腐蚀 会 严重 影响 材料 的 性 能 , 给 ADS 
的 安全 性 带 来 隐患。 故而 ,，LBE 对 容器 和 结构 钢 的 
腐蚀 问题 被 视 为 研究 ADS 系统 问题 的 焦点 之 一 . BE 
别 是 在 高 温 下 ， 这 种 液态 材料 对 铁 制 容器 有 极 强 的 
腐蚀 性 69, 虽然 有 一 些 实验 研究 了 LBE 的 腐蚀 问 
题 6- 引 , 但 是 关于 腐蚀 的 机 理 尚 未 有 透彻 的 理解 . 就 
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目前 来 说 , 腐蚀 问题 是 ADS 系统 设计 与 运行 的 一 个 
至 关 重 要 和 充满 挑战 的 问题 ,关乎 到 整个 系统 的 安 
全 性 与 可 靠 性 。 因 此 需要 从 微观 的 角度 理解 原子 间 
的 扩散 过 程 ， 从 根本 上 掌握 腐蚀 产生 的 原因 以 保证 
ADS 系统 安全 可 靠 稳定 的 运行 . 虽然 实验 结果 可 靠 ， 
但 是 其 花费 较 高 且 在 ADS 系统 中 进行 较为 困难 . 
此 通过 分 子 动力 学 模拟 来 研究 此 问题 成 为 一 个 很 好 
的 选择 . Artoto, Wang 等 人 19 也 用 分 子 动力 学 的 
方法 研究 了 铁 在 液态 馈 中 的 扩散 系数 . Manulana 等 
AUI 用 分 子 动力 学 的 方法 研究 了 Pb. Bi 原子 渗透 
到 铁 中 的 深度 。 本 文 把 腐蚀 问题 看 作为 铁 自 身 结构 
被 破坏 且 铁 原子 与 LBE 原子 相互 扩散 的 过 程 . 为 了 
揭示 铁 在 LBE 中 的 腐蚀 机 理 和 利用 在 铁 中 加 入 金属 
4 (Cr). 98 (Ni) 来 减轻 该 腐蚀 ,本 文通 过 计算 铁 及 
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其 他 元 素 原子 在 LBE 中 的 扩散 系数 , 并 得 到 了 铁 在 
LBE 中 腐蚀 过 程 中 的 微观 结构 图 并 进行 了 深入 的 分 
pr. 
1 理论 背景 

分 子 动力 学 模拟 是 一 种 用 来 计算 一 个 经 典 体 系 
的 平衡 和 传递 的 方法 ,在 特定 的 模拟 条 件 下 (温度 ， 
EJ, 时 间 ， 原子 数 等 ), 追踪 和 计算 出 原子 的 轨迹 
从 而 得 到 系统 各 种 性 质 。 在 原子 系统 中 ， 原 子 之 间 
的 相互 作用 力 是 基于 牛顿 运动 方程 : 

2 N 
M X Fj (1) 


dt? 
j=1,j#1 


AP, Fy 表示 原子 i 与 原子 ;之 间 的 相互 作用 力 , M 
是 原子 i 的 质量 , ri 表示 原子 i 的 位 移 , N 表示 原子 
总 数 ，dt 表示 模拟 时 间 步 长 。 

在 分 子 动力 学 模拟 中 ， 原 子 之 间 相 互 作用 的 势 
函数 的 选取 尤为 重要 ， 在 模拟 中 ， 原 子 之 间 采 用 
Lennard-Jones (LJ) n — m S8 gg: 


Mee 9T 
公式 2) 中 o 和 = JE LI 势 函数 参数 , 分别 表示 尺寸 
参数 、 势 阱 深度 . 

在 本 文中 ,我 们 取 n= 12, m = 6, 则 LJ AR 
数 的 形式 为 : 


va-e|("- Gy] ^ 


对 于 同 种 元 素 的 原子 (Fe-Fe, Pb-Pb, Cr-Cr, Ni-Ni) 
之 间 的 LJ 势 参数 ， 我 们 选取 Zhen 和 Davies?! 给 
出 的 值 . 饮 元 素 原 子 (Bi-Bi) 之 间 的 LJ 势 参 数 源 于 
Lemmon 等 人 [3 的 研究 . 不 同 原子 之 间 (Fe-Pb、 Fe- 
Bi, Pb-Bi, Pb-Cr, Bi-Cr, Pb-Ni, Bi-Ni) 的 势 函 数 则 
根据 Lorentz-Berthelot 混合 原则 : 


Tij = (cà + 07;)/2 (4) 
Eij = VE Ejj (5) 


本 文采 用 的 LJ 势 参数 见 表 1. 

一 般 而 言 ， 腐 蚀 是 一 个 化 学 反应 的 过 程 。 然 
m, Manly 研究 发 现 ， 高 温 时 铁 在 液态 金属 中 的 
腐蚀 , 腐蚀 过 程 中 并 未 发 现 电 子 在 原子 间 的 转移 , 非 
化 学 行为 占 了 主要 的 因素 . 因此 , 本 文 把 腐蚀 过 程 看 
作为 一 个 物理 过 程 , 而 并 非 化 学 过 程 . 本 文 主要 关注 
的 是 铁 在 LBE 中 的 扩散 过 程 , 并 应 用 以 下 两 个 公式 
计算 铁 在 LBE 中 的 扩散 系数 : 


MSD = (|F (t) —7(t= 0|) (6) 


X (6). (7) 中 MSD 表示 均 方 根 位 移 ,“()” 表 
示 系 综 平均 , 7(t) 和 7(t = 0) 分 别 表示 原子 在 时 间 t 
M t= 0 时刻 的 位 移 . 


表 1 LJ 参数 表 
Table 1 The LJ potential parameters 
元 素 e/eV o/10-10m 
Fe-Fe 0.5193 2.3193 
Pb-Pb 0.2363 3.1888 
Bi-Bi 0.0590 3.0500 
Cr-Cr 0.4989 2.3357 
Ni-Ni 0.5190 2.2808 
Fe-Pb 0.3503 3.6979 
Fe-Bi 0.1750 3.5369 
Pb-Bi 0.1181 4.8629 
Fe-Cr 0.5090 2.7086 
Fe-Ni 0.5191 2.6449 
Pb-Cr 0.3433 3.7240 
Pb-Ni 0.3502 3.6365 
Bi-Cr 0.1716 3.5619 
Bi-Ni 0.1750 3.4782 
Cr-Ni 0.5088 2.6636 


2 模拟 结果 与 讨论 分 析 


为 了 研究 铁 在 LBE 中 的 腐蚀 以 及 加 入 金属 
Cr, Ni 在 该 腐蚀 过 程 中 的 作用 ,本 文 研究 内 容 主要 
分 为 下 面 的 几 个 部 分 。 
2.1 不 同 温度 下 铁 在 LBE 中 的 扩散 系数 计算 

本 文采 用 的 模拟 盒子 为 : 中 央 建 立 17 x 17 x 17 
的 品格 常数 为 2.8682 x 107! nm 的 BCC 结构 的 铁 
原子 , 其 个 数 为 10745; 铁 原子 周围 为 LBE 原子 , 其 
中 铅 原子 个 数 为 21568(w(Pb)=0.45), 馈 原 子 个 数 为 
26137(w(Bi)=0.55), 铁 、 铅 和 锐 元 素 的 相对 原子 质量 
分 别 为 55.845, 207.2, 208. 98; 模拟 盒子 的 体积 大 小 
为 12.26 x 12.26 x 12.26 nm3,， 其 结构 如 图 1 所 示 . 对 
于 ADS 而 言 ， 其 工作 温度 范围 为 773~873 K, 我 们 
不 仅 要 研究 正常 工作 温度 下 情况 ， 而 且 需 要 考虑 极 
端 高 温 的 情况 , 因此 , 在 计算 铁 在 LBE 中 扩散 系数 
时 , 选取 温度 为 : 773 K, 873 K, 973 K, 1023 K, 1073 
K, 1173 K, 1273 K, 本 文 使 用 LAMMPS 软件 进行 模 
拟 时 所 采用 的 参数 : NVT 系 综 ， 时 间 步 长 为 0.0001 
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ps, SERERE [RIS 10 ps, 平衡 时 间 为 80 ps. 模拟 完成 
后 , 根据 公式 (6). (7) 计算 出 铁 在 LBE 中 的 扩散 系 
数 . 


NN NV M 


LBE 


图 1 Fe-LBE 系统 微观 结构 图 (Fe, Pb. Bi 原子 颜色 分 别 为 淡 
黄色 、 红色 、 蓝 色 ) 

Fig. 1 The simulation box of Fe-LBE system (The colors of 

Fe, Pb and Bi atoms are yellow, grey and indigo, respectively) 


利用 铁 在 LBE 中 的 MSD( 图 2) 和 式 (6). (7) 
计算 出 了 1023 K 时 铁 在 LBE 中 的 扩散 系数 为 
2.24 x 107? m?/s; 该 值 与 Banerje05 给 出 的 实验 值 
2.27 土 0.11 x 107? m?/s 十 分 接近 , 这 说 明了 本 文 模 
拟 结 果 的 可 靠 性 . 紧 接 着 ， 本 文 计算 了 在 不 同 温度 
下 (773 K, 873 K, 973 K, 1073 K, 1173 K, 1273 K) 
ERTE LBE 中 的 扩散 系数 , 结果 如 图 3 所 示 . 到 目前 
为 止 , 仅 Banerje 得 出 温度 为 1023 K 时 铁 在 LBE 中 
的 扩散 系数 实验 值 ， 因 此 本 文 模拟 所 得 到 的 值 可 以 
提供 有 效 的 参考 数据 . 


MSD/10-"m) 


0 20 40 60 80 
tips 


图 2 1023 K 时 铁 在 LBE 中 的 MSD 曲线 图 
Fig. 2 The MSD curve of the iron in the LBE at 1023K 
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图 3 不 同 温度 下 铁 在 LBE 中 的 扩散 系数 
Fig. 3 The diffusion coefficient of the iron in the LBE at 


different temperatures 


2.2 铁 中 挫 人 不 同比 例 Cr, Ni 

基于 上 述 模拟 的 准确 性 ， 本 文 研究 了 在 铁 中 加 
入 不 同比 例 Cr, Ni 的 情况 ， 此 时 模拟 温度 ,模拟 参 
数 都 与 2.1 节 相 同 . 本 文 在 铁 中 掺 入 金属 Cr, Ni 质 
量 分 数 如 下 : 

1) 当 铁 中 掺 入 质量 分 数 分 别 为 w(Cr) = 
596, w(Cr) = 1096, w(Cr) = 1596 的 铬 时 ， 计 算 结 
果 如 图 4 Brzn. 


一 一 无 杂质 

- -- - a(Gr)-596 
Ho o(Gr)-1096 
--- (Gr)=15% 


D/07?m?/s) 


773 873 973 1073 — 173 1273 
TIK 

图 4 不 同 温度 下 铁 中 加 入 不 同比 例 Cr 时 其 在 LBE 中 的 
扩散 系数 


Fig. 4 The diffusion coefficient of the iron in the LBE after 
adding different proportion of Cr in the iron at different 


temperatures 


2) 当 铁 中 掺 入 质量 分 数 分 别 为 w(Ni) = 596, 
w(Ni) = 1096, w(Ni) = 1596 的 镍 时 ,计算 结果 如 图 
5 Br. 
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一 一 无 杂质 
--- e(Ni)-596 
| e(Ni)-1090 
---- e(Ni)-1596 


D/10?m?/s) 


E I 873 973 103 17 1273 
TK 
图 5 不 同 温度 下 铁 中 加 入 不 同比 例 Ni 时 其 在 LBE 中 的 扩散 
系数 


Fig. 5 The diffusion coefficient of the iron in the LBE after 
adding different proportion of Ni in the iron at different 
temperatures 


3) 当 铁 中 同时 掺 入 质量 分 数 分 别 为 o (Cr) = 
596. w(Ni) 2596, w(Cr) 21096, w(Ni)=10%, w(Cr)= 
1596, (Ni) 21576 的 铬 锦 时 ,计算 结果 如 图 6 所 示 . 

ERE LBE 中 的 扩散 系数 越 大 , 说 明 其 受到 的 腐 
蚀 程 度 越 重 . 图 4 中 ,， 铁 的 扩散 系数 随 着 温度 的 升 
高 而 升 高 , 这 说 明 高 温 下 铁 腐 蚀 程 度 比 低温 下 要 大 ; 
在 不 同 温度 下 , 铁 中 加 入 不 同比 例 Cr 元 素 后 , 其 在 
LBE 中 扩散 系数 都 比 纯 铁 的 小 ， 这 说 明 Cr 的 加 入 
使 得 铁 的 腐蚀 程度 得 到 了 减轻 ; 随 着 Cr 元 素 比例 的 
增加 ， 扩 散 系 数 也 逐渐 减 小 ， 这 说 明 在 加 入 质量 分 
数 为 5%~15% 的 Cr 时 ， 随 着 Cr 质量 增加 ， 铁 腐蚀 
程度 逐渐 减 小 , 但 是 Cr 比例 在 10%、15% 时 ， 对 铁 
的 腐蚀 程度 并 不 是 很 明显 . 如 图 5 所 示 ， 铁 中 加 入 
Ni 原子 后 也 与 上 述 规律 类 似 . 而 图 6 中 , 在 铁 中 同 


时 加 入 不 同比 例 Cr、Ni 时 , 该 扩散 系数 的 减 小 速度 
比 单独 加 入 Cr, Ni 时 要 快 , 尤其 是 在 高 温 区 ,这 说 
够 更 好 地 抑制 铁 的 腐蚀 . 


明 同 时 加 入 Cr、Ni 能 


--—-q(Gr-596. e(Ni)-596 
Has e(Gr)-1096.. o(Ni)-1096 
-+= e(Gr)-1596. o(Ni)-1596 


7n 873 973 1073 1173 1273 
T/K 
图 6 不 同 温度 下 铁 中 加 入 不 同比 例 Cr, Ni 时 其 在 LBE 中 的 扩 


Fig. 6 The diffusion coefficient of the iron in the LBE after 
adding different proportion of Ni and Cr in the iron at 
different temperatures 


通过 分 析 铁 在 LBE 中 的 构 型 图 (图 7), 可 以 发 
现 加 入 金属 Cr Jet, 扩散 到 铁 中 的 LBE 原子 比 未 加 
入 时 要 少 , 而 且 随 着 加 入 比例 的 增加 而 逐渐 减少 , 这 
说 明 加 入 金属 Cr 能 阻止 铁 原 子 与 LBE 原子 相互 渗 
nee he 综 上 分 析 ， 本 文 得 出 金属 

、Ni 能 有 效 保 护 铁 的 微观 机 制 为 : 金属 Cr. Ni 
mu LBE 原子 渗入 铁 材 料 中 破坏 其 
结构 . 


(c) w(Cr)=10% 
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(b) w(Cr)=5% 
图 7 1173 K 时 铁 在 加 入 不 同比 例 Cr 时 在 LBE 中 的 微观 结构 图 (Cr 原子 为 深 黄色 ) 


Fig. 7 The microstructure of the system after adding different proportion of Cr in the iron at 1173 K (The color of Cr atom is 


deep yellow) 


结 论 


本 文采 用 分 子 动力 学 方法 研究 了 铁 在 LBE 中 的 


腐蚀 问题 , 计算 得 到 的 铁 在 LBE 中 的 扩散 系数 可 以 


提供 有 效 的 参考 . 研究 发 现 , 铁 在 LBE 中 被 腐蚀 的 


微观 机 制 是 由 于 铁 原子 与 LBE 原子 相互 渗透 , 破坏 
了 铁 材料 结构 RT, BERBA Cr, Ni 可 以 有 效 减 
轻 铁 的 腐蚀 ， 而 且 随 着 加 入 比例 的 增加 腐 乌 程度 随 
之 减 小 , 其 微观 机 理 为 金属 Cr、Ni 加 入 到 铁 后 阻止 
T LBE 原子 渗入 铁 材料 中 破坏 其 结构 ， 
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